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The  final  report  for  Air  Force  Grant  No.  F-49620- 
79-C-0193  consists  of  the  annual  progress  reports  datpd 
November  21,  1979  and  August  25,  1980,  as  well  as  the 
following  report  for  1980-81: 

During  1980-81  various  numerical  experiments  were 
made  on  one-dimensional  model  inverse  problems  modelled 
by  a  method  that  is  extendable  to  higher  dimensions.  The 
underlying  problem  is  to  recover  the  speed  of  propagation 
or  the  shape  of  an  object  from  a  scattered  field. 

The  one  dimensional  problem  that  was  investigated 
was  based  on  the  wave  equation  with  potential 
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where  the  function  to  be  recovered  is  the  potential  q(x). 
The  model  initial  conditions  were  u  =  0  ,  u.  *  “26'  (x)  . 
The  model  boundary  condition  was  u  =  0  on  x  =  0  .  And 
the  "extra"  condition  which  determines  q(x)  is  either  (a) 
3u/3x  on  x  =  0  for  all  time  or  (b)  the  scattered  field 
u'vfi  (t-x)  +  R(t-x)  at  x  ■  +  «>  .  The  free  space  (q=0) 
solution  is  u  =  6 (x-t)  -  <5 (x+t)  .  The  distorted  plane 
wave  P(x,t)  is  defined  by  the  condition  that  it  is  a 
solution  of  the  equation  and  at  x  =  •  00  ,  t  -*•  +  00  , 

82  01  £8  03  2 


Approved  for  public  rcloaso  , 
distribution  unlimited. 


-2- 


it  becomes  6(t-x)  (exactly  if  q  has  compact  support) 

and  P  is  written  as  6(t-x)  +  P  .  . 

scat 

The  linearized  equation  based  on  the  assumption 

2 

that  q<<l  and  that  terms  of  ordc  r  q  may  be 
neglected,  yields  for  (a) 

«  -  +  2  lx  (2x)  -  Pscat  (0'-2x)] 
t 

Wlth  \  H  dt 
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for  (b)  . 
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for  t  >  0  ; 


The  numerical  method  is  based  on  the  we.i  1-known 
formula  that  follows  from  the  propagation  of  singularities 


q  =  2 


d  lim 
dx  t-*-x+ 


u(x,x) 


and  then  representing  u(x,x)  by  means  of  distorted 
plane  waves.  (All  tnese  have  analogous  formulas  in 
higher  dimensions . ) 

The  analogues  of  the  above  equation  ere  for  (b) 


<5  “  ~2  lx  R(2x)  +  lx  J  (Pscat(x'x_s)“Pscat("x'x_s))R(s)  ds 
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R(s)  =  +  pscat<°'~s>  + 


psoat(0'-s',ds' 


Psoat(0't-S)dt 


The  idea  is  to  use  one  of  these  formulas  as  part  of 
an  iteration.  We  guess  q  and  P  and  recompute  q  using 
the  above  formulas  for  cases  (a)  or  (b) . 

Note  that  if  the  given  data  vanishes  in  case  (b)  we 
recover  q  =  0  but  in  case  (a)  we  donot. 

We  have  been  unable  to  implement  (a)  as  planned. 

The  reason  is  that  we  must  operate  in  a  finite  region  and 
use  a  radiation  condition  at  a  finite  distance.  In  our 
case  for  example  with  q  a  simple  quadratic  with  support 
in  |x|  <  2  ,  over  “8  <_  x  8  and  0  <_  t  <_  . . .  The 

errors  produced  are  compounded  in  case  (a) .  The  potential 
could  only  be  recovexed  to  40%  of  its  original  value  with 
an  error  of  about  5%  of  the  maximum  potential.  However, 
the  potentials  were  well  beyond  the  range  of  the  linear 
theory. 

In  case  (b)  even  with  a  coarse  mesh  the  results  are 
surprisingly  worse. 

The  results  are  being  written  up  as  a  technical 
report  to  be  issued  at  the  Courant  Institute. 
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